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PREFACE 


The  Alcatraz  disposal  modeling  work  described  in  this  report  was  per¬ 
formed  for  the  US  Army  Engineer  District,  San  Francisco  (SPN),  Planning/ 
Engineering  Division,  as  part  of  the  District's  Dredged  Material  Disposal 
Management  Program. 

The  numerical  simulations  portion  of  the  study  was  conducted  in  the 
Hydraulics  Laboratory  (HL)  of  the  US  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  (WES)  during  the  period  October  1986  to  March  1988  under  the  general 
supervision  of  Messrs.  Frank  A.  Herrmann,  Chief,  HL;  Richard  A.  Sager,  Assis-, 
tant  Chief,  HL;  William  H.  McAnally,  Jr.,  Chief,  Estuaries  Division;  and 
William  D.  Martin,  Chief,  Estuarine  Engineering  Branch. 

The  HL  work  was  performed  by  Mmes.  Paula  G.  Kee  and  Virginia  R.  Pankow, 
Estuarine  Engineering  Branch,  under  the  technical  guidance  of 
Messrs.  Michael  J.  Trawle ,  Waterways  Division,  and  Allen  M,  Teeter,  Estuarine 
Processes  Branch,  Estuaries  Division,  former  and  present  Project  Managers. 

Dr.  Billy  H.  Johnson,  Math  Modeling  Croup,  Waterways  Division,  and 
Ms.  Barbara  P.  Donnell,  Estuarine  Simulation  Branch,  Estuaries  Division,  pro¬ 
vided  advice  and  assistance  with  the  DIFID  and  TABS-2  models.  Additional  mod¬ 
eling  support  was  provided  by  Mr.  Julian  M.  Savage,  Estuarine  Processes 
Branch.  Appendix  A  was  written  by  the  HL  staff,  and  Appendix  B  was  written  by 
Mr.  Trawle  and  Dr.  Johnson.  This  report  was  prepared  by  Ms.  Pankow  and  edited 
by  Mrs.  Marsha  C.  Gay,  Information  Technology  Laboratory,  WFS. 

Mr.  Brian  Walls  was  SPN  liaison  during  the  study,  and  Mr.  Thomas 
Wakeman,  Model  Director,  SPN,  supervised  the  work  performed  at  the  San 
Francisco  Bay-Delta  Model.  Mr.  Lester  Tong,  SPN,  is  Program  Manager  of  the 
Dredged  Material  Disposal  Management  Program. 

COL  Dwayne  G.  Lee,  EN,  is  the  Commander  and  Director  of  WES. 

Dr.  Robert  W.  Whalin  is  the  Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

acres 

cubic  feet 
cubic  yards 
feet 

miles  (US  nautical) 

miles  (US  statute) 

pounds  (force)- 

second  per  square  foot 

square  feet 


_ By. _ 

4,046.8/3 

0.02831685 

0.7645549 

0.3048 

1.852 

1.609347 

47.88026 

0.09290304 


_ To  Obtain 

square  metres 
cubic  metres 
cubic  metres 
metres 
kilometres 
kilometres 
pascal s-second 

square  metres 
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SAN  FRANCISCO  BAY:  MODELING  SYSTEM  FOR  DREDGED  MATERIAL 


DISPOSAL  AND  HYDRAULIC  TRANSPORT 


PART  I:  INTRODUCTION 


Background 


1.  The  Alcatraz  disposal  site  near  Alcatraz  Island  is  a  naturally  deep 
open-water  site  in  San  Francisco  Bay  (Figure  1).  This  has  been  used  as  a 
dispersive  site  for  dredged  material  since  the  turn  of  the  century.  As  such, 
material  disposed  there  is  intended  to  be  transported  out  of  the  site  toward 
the  Golden  Gate  Bridge  and  out  to  sea  by  the  strong  tidal  currents. 

2.  In  1982  personnel  of  the  US  Army  Engineer  District  (USAED),  San 
Francisco,  discovered  that  the  dredged  materials  being  disposed  at  the  Alca¬ 
traz  site  were  not  being  dispersed  at  the  rate  that  they  were  being  introduced 
and  that  a  permanent  mound  was  forming.  This  mound  was  a  hazard  to  navigation 
and  indicated  that  the  capacity  of  the  disposal  site  was  being  taxed.  There 
was  also  concern  that  more  expensive  alternate  disposal  sites  would  have  to  be 
developed  in  the  future. 

I.  In  ’•esponse  to  the  problem,  USAED.  San  Francisco,  formed  the  Dredged 
Material  Disposal  Management  Program  (DMP) ,  which  was  tasked  to  (a)  determine 
the  fate  of  disposed  material  at  the  Alcatraz  disposal  site,  and  (b)  develop 
modeling  tools  to  assess  disposal  of  dredged  materials  at  various  other  loca¬ 
tions  in  San  Fr«nci.siu  D„,  .  A  request  was  made  to  the  US  Army  F.neineer  Water- 
wavs  Experiment  Station  for  assistance  witii  these  goals. 

Objective 


4.  The  overall  objective  of  the  disposal  modeling  study  is  to  assist 
USAED,  San  Francisco,  to  meet  the  goals  of  the  DMP  by  the  development  of  a 
modeling  tool  that  can  define  the  fate  of  dredged  material  disposed  at  the 
Alcatraz  disposal  site.  The  purpose  of  this  report,  is  to  describe  and  summar¬ 
ize  the  activities  involved  in  the  development  and  validation  of  the  models 
used  in  this  study  and  to  demonstrate  the  capability  of  STUDH  (Sediment 


Transport  in  Unsteady  Two-Dimensional  Flows,  Horizontal  Plane),  the  sediment 
transport  model. 

Approach 

5.  The  study  objective  was  met  through  the  use  of  both  physical  and 
numerical  modeling  techniques.  This  San  Francisco  Pay  hybrid  modeling  ap¬ 
proach  consisted  of  the  following  interrelated  parts: 

a.  The  creation  of  a  large-scale  or  comprehensive  numerical  model 
of  the  San  Francisco  Bay  to  simulate  the  hydrodynamic  character¬ 
istics  of  the  entire  bay  system.  This  numerical  model  was  used 
to  define  tidal  hydrodynamics  and  overall  circulation  patterns. 

b.  The  use  of  the  San  Francisco  Bay-Delta  Model  (SFBM)  to  generate 
water-surface  and  current  velocity  data  for  the  numerical  model 
boundary  conditions  and  the  verification  data  sets. 

c.  The  development  of  a  higher  resolution  inset  numerical  model  of 
the  Alcatraz  disposal  area.  This  model  was  driven  by  the  com¬ 
prehensive  model  and  used  to  establish  water  level  and  current 
velocity  boundarv  conditions  for  use  in  the  disposal  model. 

d.  The  use  of  the  numerical  model  D'FID  ^Disposal  From  Instanta¬ 
neous  Dump)  to  simulate  barge  or  hopper  disposal  of  the  dredged 
material.  The  model  predicts  the  short-term  fate  (less  than 

1  hr)  of  disposed  material  ,  uding  the  initial  deposition 

pattern . 

e.  The  use  of  the  sediment  transput  model  STDDH  to  simulate  the 
fate  of  the  disposed  materia]  using  suspended  sediment  concen¬ 
trations  produced  by  the  dispesu J  model. 


PART  II:  DESCRIPTION  OF  THE  MODELS 


Modeling  Overview  and  Procedures 


6.  To  calculate  the  shor.--  and  long-term  transport  of  disposed  sediment 
materials,  hydrodynamic  data  are  needed  at  many  locations  in  the  flow. 
Therefore,  the  first  step  in  sediment  modeling  is  the  development  of  the  nec¬ 
essary  dynamic  current  field.  Existing  disposal  models  lack  erosion  and  de¬ 
position  features  necessary  for  long-term  predictions.  Long-term  sediment 
transport  models  lack  disposal  prediction  capability.  Therefore,  a  combina¬ 
tion  of  models  including  a  disposal  model,  a  sediment  transport  model,  and  a 
hydrodynamic  model  were  merged  in  this  study.  The  models  selected  were  DTFID, 
STl'DH,  and  RMA-2V  ;'Two-Dimensi onal  Model  for  Open-Channel  Flows).  The  latter 
two  models  are  included  in  the  TABS-2  modeling  system  supported  by  the  US  Army 
Corps  of  Engineers. 

7.  In  addition  to  the  family  of  numerical  models,  the  physical  model  of 
San  Francisco  Bay  was  used  to  generate  svnoptic  water-surface  elevation  and 
velocity  data.  These  data  would  provide  the  necessary  boundary  condition  in¬ 
formation  and  the  verification  data  sets  for  the  numerical  models. 

The  Ban  Francisco  Bav-Pelta  Model 


8.  The  SFRM  is  a  fixed-bed,  distorted  scale,  tidal  hydraulic  model. 

This  physical  hydraulic  model  provides  a  means  of  reproducing  on  a  manageable 
scale  the  hvdrodvnamic  phenomena  that  occur  throughout  the  large  and  complex 
estuarine  system.  The  establishment  of  the  similarity  between  model  and  pro¬ 
totype  commonly  makes  use  ot  laws  of  similitude  expressed  hv  a  set  of  dimen¬ 
sionless  numbers.  These  include  the  Fronde,  Reynolds,  Weber,  and  C.auehv 
numbers,  which  are  the  ratios  of  inertial  forces  to  the  forces  of  gr;ri'.  , 
viscosity,  surface  tension,  and  elasticity,  respectively.  Perfect  similitude 
would  require  that  each  of  these  numbers  be  .lie  same  for  model  and  prototype. 
Physical  hydraulic  models  cannot  be  designed  to  provide  similitude  with  re¬ 
spect  to  all  forces,  but  are  designed  to  simulate  the  dominant  forces  affect¬ 
ing  the  conditions  being  modeled.  The  scale  is  selected  so  that  the  remaining 
forces  are  negligible  or  do  not  cause  fundamental  dissimilarities  in  the  pro¬ 
cesses  of  interest. 
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9.  In  the  San  Francisco  Bay  complex,  the  depth,  surface  slope,  and 
other  features  of  tidal  flow  are  controlled  by  the  joint  effects  of  inertial 
and  gravitational  forces.  The  major  hydraulic  quantities  vary  according  to 
the  Froude  number.  For  the  SFBM,  different  model  scales  were  chosen  for  hor¬ 
izontal  and  vertical  lengths  to  reduce  construction  and  operation  time  and 
cost  as  well  as  keep  faithful  reproduction  of  flows  and  permit  measurements  of 
the  desired  parameters  with  satisfactory  accuracy.  The  effect  of  distortion 
on  boundary  roughness  was  reproduced  during  verification  tests  by  the  use  of 
metal  strips  embedded  in  the  model.  These  strips  were  adjusted  to  change  the 
frictional  resistance  to  obtain  reproduction  of  the  bay's  tides,  current,  and 
salinities.  The  scaling  ratios  of  model  to  prototype  are  as  follows: 


Characteristic _ 

Horizontal  length 
Vertical  length 
Time 

Velocity  (horizontal) 

Discharge 

Salinity 


Scale  Relations 
Model : Prototype 

1:1,000 

1:100 

1:100 

1  :  10 

1 : 1,000,000 
1  :  1 


10.  The  SFBM  occupies  an  area  of  about  1  acre,*  and  is  completely  en¬ 
closed  in  a  1 28 , 500-sq-f t  shelter  to  protect  it  from  the  weather  and  thus  per¬ 
mit  uninterrupted  operation.  The  limits  of  the  model,  shown  in  Figure  2, 
encompass  a  portion  of  the  Pacific  Ocean  extending  17  miles  west  of  the  Golden 
Gate,  San  Francisco  Bay,  San  Pablo  Bay,  Suisun  Bay,  and  all  of  the  Sacramento- 
San  Joaquin  Delta  east  of  Suisun  Bay  to  the  cities  of  Sacramento  on  the  north, 
Stockton  on  the  east,  and  Tracy  to  the  south. 

11.  The  model  is  equipped  with  the  necessary  appurtenances  including 
freshwater  and  saltwater  supply  sumps,  primary  saltwater  inflow  and  mixing 
pumps,  computer-controlled  tide  generator  capable  of  reproducing  repetitive 
ocean  tides  by  regulating  two  motorized  gate  valves,  precise  flowmeters  to 
regulate  major  river  inflows  to  the  delta  plus  other  inflow  and  withdrawal 
facilities,  major  pumping  plant  exports,  and  electrically  operated  one-way 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  3. 
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head  gates.  The  model  represents  the  state  of  the  art  in  model  operation, 
real-time  data  acquisition  using  miniature  sensors,  data  analysis,  anc  infor¬ 
mation  storage/retrieval  capability. 

The  TABS-2  Modeling  System 

12.  The  TABS-2  is  a  modular  system  composed  of  distinct  computer  pro¬ 
grams  linked  together  by  preprocessors  and  postprocessors.  Each  of  the  major 
computer  programs  solves  a  particular  type  of  problem:  hydrodynamics 
(RMA-2V) ,  sediment  transport  (STUDH) ,  or  water  quality  (RMA-4) .  These  pro¬ 
grams  employ  the  finite  element  method  to  solve  the  governing  equations.  A 
brief  description  of  RMA-2V  and  STUDH  appears  in  Appendix  A. 

13.  The  major  programs  from  the  TABS-2  system  used  in  this  study  are 
RMA-2V  and  STUDH.  Along  with  a  disposal  model,  DIFID,  the  three  programs 
interact  as  shown  in  the  following  diagram: 

COMPREHENSIVE  INSET  DISPOSAL 

RMA-2V  RMA-2V  DIFID 

hydrodynamic - ►hydrodynamic  ►hvdrodvnamic 

I 

INSET 

STUDH  DTFTD 

sediment  sediment 

transport  «  initial  concentration 

14.  A  summary  of  program  codes  used  in  this  modeling  effort  and  a 

flowchart  of  the  code  usage  are  presented  in  Table  1  and  Figure  3,  respec¬ 
tively.  The  procedures  of  model  development  and  execution  were  followed  as 
described  in  the  TABS-2  user's  manual  (Thomas  and  McAnally  1985). 

Hydrodynamic  Model 

15.  RMA-2V  is  a  finite  element  solution  of  the  Reynolds  form  of  the 
Navier-Stokes  equations  for  turbulent  flows.  Friction  is  calculated  with 
Manning's  equation,  and  eddy  viscosity  coefficients  are  used  to  define  turbu¬ 
lent  exchange  characteristics.  A  velocity  form  of  the  basic  equation  is  used 
with  side  boundaries  treated  as  either  slip  (parallel  flow)  or  static  (zero 
flow).  The  model  recognizes  computationally  wet  or  dry  elements  and  corrects 
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Program 

GFGEN 

EDGRG 

CODE24 

RMA-2V 

JOBSTRM 

VPLOT 

POSTHYD 

DIFID 

ENGMET 

STUDH 

SEDGRAF 


Table  1 

Summary  of  Program  Codes 


_ Description _ _ _ _ 

Create  geometric  data  file  for  use  by  models 

Edit  computational  network  interactively 

Process  physical  model  data  files  for  input  to  RMA— 2V  for  use  as 
time-dependent  boundary  conditions 

Compute  two-dimensional  vertically  averaged  flows  and  water 
levels 

Compute  boundary  conditions  for  an  inset  computational  network 
from  results  of  a  larger  network  model  run 

Produce  vector  plots  of  RMA-2V  velocity  results 

Analyze  results  of  an  RMA-2V  run  and  plot  time-histories 

Compute  suspended  sediment  concentrations  and  initial  deposition 
on  the  bottom  of  material  released  into  the  water  column  from  a 
dredged  material  disposal  operation 

Translate  RMA-2V  output  from  non-SI  (English)  units  to  SI 
(metric)  units  for  use  by  STUDH 

Compute  two-dimensional  vertically  averaged  transport,  erosion, 
and  deposition  of  sediments 

Produce  factor  map  of  STUDH  results 
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GFGEN 


the  mesh  accordingly.  Boundary  conditions  may  be  water-surface  elevations, 
velocities,  or  discharges  and  may  occur  inside  the  mesh  as  well  as  along  the 
outer  boundaries. 

Comprehensive  mesh 

16.  Model  limits.  The  comprehensive  or  global  model  was  established  as 
extending  from  7  nautical  miles  west  of  the  mouth  of  San  Francisco  Bay  to 
Carquinez  Strait  between  San  Pablo  and  Suisun  bays.  A  mesh  was  developed  con¬ 
sisting  of  1,580  quadrilateral  elements  with  5,000  nodes.  The  bathymetric 
contours  of  the  bay  were  major  factors  in  the  location,  shape,  and  size  of  the 
elements.  The  boundaries  were  represented  by  curved  element  sides,  ar^  the 
model  included  the  five  islands  of  Alcatraz,  Angel,  Brooks,  and  Treasure/Yerba 
Buena  and  lower  San  Pablo  Bay. 

17.  Bathymetry.  Bathymetric  data  input  to  the  model  were  derived  from 
the  following  1 :40, 000-scale  National  Ocean  Survey  charts: 


Chart  No. 

Location 

Date 

18651 

San  Francisco  Bay — Southern  Part 

1978 

18649 

Entrance  to  San  Francisco  Bay 

1985 

18654 

San  Pablo  Bay 

1983 

18656 

Suisun  Bay 

1983 

18.  An  x-coordinate ,  y-coordinate ,  and  z-value  (bed  elevation)  were 
assigned  to  each  of  the  5,000  nodes  during  the  digitizing  process  and  stored 
in  a  computer  file  for  further  use  and  modification.  The  vertical  datum  was 
mean  lower  low  water  (ml lw) .  Element  types  were  established,  and  eddy  viscos¬ 
ity  (turbulent  exchange  coefficient)  and  Manning's  n  (bottom  friction)  val¬ 
ues  were  assigned  to  each  element  type. 

19.  The  network  generation  program  GFGEN  was  used  to  create  a  geometry 
and  finite  element  computational  mesh  file  for  input  to  the  TABS-2  modeling 
system.  A  utility  program,  EDGRG,  was  used  to  interactively  edit  the  computa¬ 
tional  mesh.  With  EDGRG,  curved  boundaries  can  be  fitted  and  nodes  and  ele¬ 
ments  of  the  mesh  can  be  moved,  added,  or  deleted.  These  changes  can  be  saved 
as  a  permanent  file. 

20.  Boundary  conditions.  If  no  boundary  condition  is  specified  at  a 
node,  the  program  computes  velocity  x-  and  y-components  and  the  water  depth  at 
that  node.  For  the  nodes  along  solid  boundaries  (land/water  interface),  slip 
flows  parallel  to  the  boundary  were  specified.  A  head  specification  (water- 
surface  elevation)  was  assigned  to  the  Pacific  Ocean  boundary  nodes.  These 
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nodes  had  a  water-surface  elevation  specified  at  each  time-step.  The  fresh¬ 
water  inflow  boundaries  at  Benicia,  Mare  Island,  and  South  Bay  had  velocity 
specifications  (Figure  4) . 

21.  The  input  data  for  the  boundary  conditions  of  the  numerical  model 
were  supplied  by  the  SFBM  (USAED,  San  Francisco,  in  preparation).  Four 
conditions  were  simulated  in  the  physical  model: 

a.  Nineteen-year  mean  ocean  tide  with  low  freshwater  delta  outflow 
(4,400  cfs). 

b.  Nineteen-year  mean  ocean  tide  with  high  freshwater  delta  out¬ 
flow  (40,000  cfs). 

£.  Spring  ocean  tide  with  low  freshwater  delta  outflow. 

d.  Spring  ocean  tide  with  high  freshwater  delta  outflow. 

Synoptic  data  were  taken  at  10  tide  and  36  velocity  stations  under  each  of  the 
four  test  conditions.  These  data  were  used  to  provide  boundary  conditions  and 
verify  the  numerical  hydrodynamic  model. 

22.  Conversion  of  the  physical  model  data  as  input  to  the  RMA-2V 
hydrodynamic  model  for  use  as  time-dependent  boundary  conditions  was  performed 
by  the  utility  code  CODE24.  CODE24  reads  files  containing  physical  model  ve¬ 
locities  and  water-surface  elevations  and  produces  a  boundary  control  update 
file  for  RMA-2V.  This  file  contains  time-varying  water-surface  elevations  and 
x-  and  y-direction  velocity  components. 

23.  Time-step .  The  time-step  used  for  the  numerical  simulations  was 
the  same  as  that  used  for  data  collection  on  the  physical  model.  The  24.84-hr 
lunar  day  was  divided  into  40  equal  parts  producing  a  time-step  of  37.26  min¬ 
utes.  The  same  time  increment  was  used  to  make  direct  physical  model/ 
numerical  model  data  comparisons. 

24.  RMA-2V  output ■  The  output  files  from  an  RMA-2V  run  consist  of  a 
water-surface  elevation  and  velocity  file  in  binary  format,  a  binary  hot-start 
file  (which  can  be  used  to  start  a  continuation  run) ,  and  a  line  printer  ASCII 
output.  No  graphics  are  produced  by  RMA-2V.  Instead  the  output  file  is  ready 
for  use  by  utility  programs,  which  prepare  velocity  vector  plots,  time-history 
plots,  contours,  or  other  plotted  outputs. 

Inset  mesh 

25.  An  area  containing  Alcatraz  Island  and  the  Alcatraz  dredged  mate¬ 
rial  disposal  site  was  chosen  to  become  the  higher  resolution  "inset"  mesh 
(Figure  5).  The  disposal  area,  which  consisted  of  one  element  in  the  global 
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mesh,  was  changed  to  25  small  higher  resolution  elements  in  the  inset  mesh. 
All  elements  in  the  lower  Central  Bay  and  Golden  Gate  areas  were  reduced  in 
size  to  increase  resolution  in  the  area  affected  by  the  disposal  operations. 
The  digitized  mesh  contained  x-  and  y-coordinates  and  bed  elevations  for  584 
elements  and  1,802  nodes.  An  interface  program,  JOBSTRM,  was  used  to  create 
an  inset  mesh  boundary  condition  file  for  RMA-2V  using  output  files  from  pre¬ 
vious  runs  of  the  program  on  the  global  mesh.  Water-surface  elevations  were 
specified  at  the  ocean  boundary,  and  velocity  specifications  were  assigned  to 
the  remaining  water  boundaries  (Figure  6) . 

Disposal  Model 


26.  The  model  used  for  the  Alcatraz  disposal  site  simulation  was  DIFID. 
This  numerical  model  accounts  for  the  physical  processes  that  determine  the 
short-term  fate  of  dredged  material  disposed  at  an  open-water  site.  Disposal 
phases  include  convective  descent,  dynamic  collapse  and  spreading,  and  long¬ 
term  diffusion.  The  long-term  diffusion  phase  in  DIFID  assumes  a  uniform  cur¬ 
rent  field,  which  is  unrealistic  for  the  Alcatraz  disposal  site.  Therefore, 
long-term  diffusion  calculations  were  performed  using  the  STUDH  sediment 
transport  model.  Output  from  DIFID  includes  an  estimate  of  suspended  sediment 
concentration  and  the  initial  deposition  of  material  on  the  bottom.  An  ex¬ 
planation  of  the  theoretical  basis  of  DIFID  (Trawle  and  Johnson  1986a)  can  be 
found  in  Appendix  B.  Grid  development  and  model  execution  were  performed 
according  to  Johnson  (in  preparation). 

Model  limits 

27.  The  DIFID  model  requires  a  grid  of  uniformly  sized  square  cells. 
Therefore,  a  grid  of  17  by  39  cells,  each  cell  400  ft  by  400  ft,  was  used  to 
represent  the  area  identified  in  Figure  7.  The  input  file  contained  a  depth 
for  each  grid  point  and  information  pertaining  to  the  disposal  site  (ambient 
current  velocities,  vertical  mixing  coefficients,  and  vertical  water  density 
gradient),  the  character  of  the  dredged  material  (bulk  density,  sediment  com¬ 
position,  and  particle  velocities)  ,  and  the  disposal  operation  (barge  position 
on  grid,  vessel  dimensions  and  velocity,  and  volume  of  material  to  be 
disposed) . 

Hydrodynamic  input 

28.  Using  the  binary  output  files  from  the  inset  mesh  RMA-2V 
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Figure  6.  Alcatraz  inset  mesh  with  boundary  condition  nodes 
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simulations,  the  TABS-2  program  VPLOT  was  used  to  create  a  special  velocity 
output  file  that  was  interpolated  onto  a  rectangular  grid.  This  rectangular 
grid  velocity  file,  which  contained  x-  and  y-location  coordinates  and  x-  and 
y-velocity  components,  was  used  to  establish  the  velocity  field  for  the  dis¬ 
posal  model.  The  DIFID  model  was  run  under  conditions  of  maximum  flood  and 
ebb  velocities. 

DIFID  output 

29.  The  user-specified  output  from  these  simulations  was  in  the  form  of 
a  computer  printout  that  listed  water  column  suspended  sediment  concentrations 
and  bottom  accumulation  for  each  solid  fraction  identified  in  the  input  file, 
at  the  designated  depths  and  times  following  the  disposal  operation. 

Sediment  Transport  Model 


30.  The  sediment  transport  model  in  the  TABS-2  system  is  STUDH.  This 
model  solves  the  convection-diffusion  transport  equation  with  bed  source 
terms.  These  terms  are  structured  for  either  sand  or  cohesive  sediments.  The 
Ackers-White  (1973)  procedure  is  used  to  calculate  a  sediment  transport  poten¬ 
tial  for  the  sands  from  which  the  actual  transport  is  calculated  based  on  the 
availability.  Clay  erosion  is  treated  as  described  in  Ariathurni,  MacArthur, 
and  Krone  (1977),  and  the  deposition  of  clay  uses  Krone's  (1962)  equations. 
Deposited  material  forms  layers,  and  bookkeeping  within  the  STUDH  code  allows 
up  to  10  layers  at  each  node  for  maintaining  separate  material  types,  depos¬ 
ited  thickness,  density,  and  age.  The  code  uses  the  same  mesh  as  RMA-2V. 

This  model  was  used  to  predict  suspended  sediment  transport  resulting  from 
sediment  introduction  into  the  water  column  during  the  disposal  operation. 
Input 

31.  Input  information  necessary  for  STUDH  includes  the  geometry  file 
from  GFGEN ,  velocity  and  water-surface  files  from  RMA-2V,  and  the  run-control 
file  containing  suspended  sediment  characteristics  and  concentration  (results 
from  the  disposal  operation  plus  any  background  concentrations  that  might  be 
present).  The  inset  mesh  geometry  was  used,  and  with  the  use  of  the  utility 
program  ENGMET,  the  GFGEN  and  RMA-2V  files  were  converted  from  non-ST  to 

SI  (metric)  units  as  required  fot  STUDH. 

STUDH  output 

32.  The  binary  output  file  from  STUDH  consists  of  sediment 
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concentration,  velocity,  water  depth,  and  bed  changes  at  designated  nodes.  No 
graphics  are  produced  by  STUDH;  however  the  output  file  is  ready  for  use  by 
utility  programs  which  prepare  plots  of  sediment  transport  and  new  bottom 
contours . 

Model  Integration 

33.  As  can  be  seen  from  the  preceding  discussion,  utility  programs 
within  the  TABS-2  system  enable  output  files  from  one  program  to  become  input 
to  other  programs  in  the  system.  Program  output  is  in  a  format  that  is  usable 
by  a  set  of  plot  programs  to  display  the  data  in  several  formats.  The  same 
mesh  can  be  used  by  the  hydrodynamic  and  the  sedimentation  transport  codes, 
saving  time  and  confusion.  Utility  programs  are  available  to  edit  and  change 
the  mesh.  This  is  economical  for  small  changes;  however,  in  this  study,  it 
was  easier  and  less  time-consuming  to  redigitize  and  renumber  the  inset  mesh. 
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PART  III:  VERIFICATION  AND  APPLICATION  OF  THE  MODELS 


Physical  Model  Verification 


34.  The  verification  report  of  the  San  Francisco  Bay  physical  model 
(USAED,  San  Francisco,  1963)  indicated  that  the  model  could  reproduce  a  proto¬ 
type  mean  tide  (21-22  September  1956)  with  good  results.  Maximum  error  in 
elevation  in  model  reproduction  at  23  gages  throughout  the  bay  was  0.4  ft  pro¬ 
totype,  with  the  model  amplitudes  being  slightly  less  than  those  of  the  proto¬ 
type.  The  verification  procedure  indicated  that  local  model  velocities  may 
vary  up  to  20  percent  from  those  in  the  prototype,  the  model  velocities  tend¬ 
ing  to  be  generally  higher  than  those  in  the  prototype.  Additional  prototype 
to  model  comparisons  made  in  1976  (USAED,  San  Francisco,  1976)  were  similar. 

Hydrodynamic  Model  Verification 


35.  Verification  of  the  RMA-2V  model  consisted  of  several  test  runs  to 
adjust  the  boundary  conditions  and  internal  coefficients  so  the  numerical 
model  would  reproduce  water-surface  elevations  and  current  velocities  supplied 
by  the  SFBM.  Verification  of  the  model  started  with  static  leak  tests  to 
ensure  the  integrity  of  the  numerical  mesh.  Next,  steady  state  and  dynamic 
conditions  were  imposed  at  the  boundaries.  Progressive  sophistication  of 
these  tests  included  adjusting  one  parameter  at  a  time  to  determine  the  effect 
on  the  model.  This  iterative  process  of  adjustments  resulted  in  the  final 
global  mesh  with  designated  internal  coefficients  and  dynamic  boundary 
conditions. 

Test  conditions 

36.  Two  sets  of  physical  model  data  were  used  tc  verify  the  hydrody¬ 
namic  model.  First  a  19-year  mean  tide  was  generated  at  the  ocean  tide  gener¬ 
ator  to  produce  a  mixed  tide  with  extremes  of  2.3,  5.8,  0.0,  and  4.6  ft  mllw 
at  sta  8,  Presidio,  Golden  Gate.  The  delta  net  outflow  was  maintained  at 
4,400  cfs.  After  all  adjustments  were  made  and  the  RMA-2V  model  was  satisfac- 
torially  reproducing  the  mean  tide  physical  model  data,  a  second  data  set  was 
used.  These  data  were  a  simulation  of  a  28-29  December  1955  spring  tide  with 
extreme  values  of  3.0,  7.2,  -1.2,  and  4.9  ft  mllw  at  sta  8  with  the  same 
^,400-cfs  delta  net  outflow. 
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Model  coefficients 


37.  The  bay  was  divided  into  four  major  eddy  viscosity  coefficient 
zones.  Each  zone  had  four  to  five  different  Manning's  n  values  depending  on 
the  water  depth  and  the  size  of  the  element,  with  greater  roughness  (0.03)  for 
shallow  areas  and  less  roughness  (0.015)  for  deep  areas.  Eddy  viscosity  coef¬ 
ficients  were  assigned  according  to  the  turbulence  characteristics  of  the 
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region  and  ranged  from  200  to  1,000  Ib-sec/ft  . 

38.  Physical  model  data  were  put  into  the  format  compatible  with  the 
TABS-2  system  requirements.  The  velocity  boundary  specifications  were  ad¬ 
justed  in  magnitude  and  direction  when  necessary  to  create  entrance  conditions 
needed  to  produce  the  desired  results  inside  the  mesh.  Similarly,  the  physi¬ 
cal  model  tide  data  at  the  ocean  tide  sta  0  was  adjusted  1  hr  in  timing  to  re¬ 
produce  the  desired  tide  at  the  tide  stations  inside  the  global  mesh. 

39.  Velocity  output  from  RMA-2V  is  vertically  averaged.  Before  com¬ 
parisons  could  be  made,  the  physical  model  data  were  also  vertically  averaged. 
Physical  model  stations  that  reported  only  middepth  values  were  unaffected, 
and  the  data  were  used  as  reported.  For  stations  that  reported  surface  (5  ft 
below  mllw)  and  bottom  (5-10  ft  above  bed)  velocity  values,  an  average  of  the 
two  was  used.  Sta  1  at  the  Golden  Cate  reported  data  from  the  surface  and 
middepth,  and  a  value  equal  to  85  percent  of  the  surface  velocity  was  used. 
Global  mesh  verification 

40.  The  numerica1  simulation  was  30  hr  with  the  data  from  the  last 

25  hr  plotted.  The  lead-in  5  hr  was  spin-up  time  and  was  used  as  a  hot-start 
input  file  for  the  final  run.  Ten  tide  and  thirteen  velocity  stations  were 
chosen  for  RMA-2V  and  physical  model  comparisons  (Figure  .  Plates  1-7  are 
tide  and  velocity  station  comparison  plots  for  the  mean  tide,  4,400-cfs  net 
delta  outflow  simulation;  and  Plates  8-14  are  the  comparisons  for  the  spring 
tide,  4,400-cfs  net  delta  outflow  conditions. 

41.  Tide  comparisons  for  mean  tide  conditions  were  good.  Phase  differ¬ 
ences  were  less  than  one  37-min  time-step,  and  amplitude  values  compared  fa¬ 
vorably.  The  tide  planes  (means)  of  RMA-2V  values  were  0.10  ft  and  0.20  ft 
lower  at  the  South  Bay  sta  2  and  Carquinez  Strait  sta  20,  respectively.  All 
interior  tide  stations  had  good  tide  plane  agreement. 

42.  Comparisons  of  EMA-2V  versus  physical  model  velocity  values  were 
reasonable.  RMA-2V  velocity  amplitudes  tended  to  be  lower  than  the  physical 
model  at  most  station.  Sta  1-4B  in  both  the  numerical  and  phys  cal  models 
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Figure  8, 


San  Francisco  Bay  area  with  physical  model  limits  and 
station  locations 


displayed  the  same  unusual  flattened  velocity  characteristics. 

A3.  Using  the  same  geometry,  eddy  viscosity  coefficients.  Manning's  n 
roughness  values,  and  net  delta  outflow,  a  model  simulation  of  the  spring  tide 
was  run.  The  same  tide  and  velocity  stations  were  examined  for  comparisons. 
All  RMA-2V  tide  and  velocity  stations  compared  favorably  with  the  physical 
model  data  and  displayed  the  same  characteristics  as  seen  in  the  mean  tide 
data  comparisons. 

Inset  mesh  verification 

AA.  The  turbulence  and  roughness  coefficients  used  in  the  inset  mesh 
were  adjusted  to  be  compatible  with  the  smaller  size  of  the  new  elements. 

Once  this  was  done,  a  simulation  of  the  mean  tide,  A,A00-cfs  net  delta  outflow 
condition  was  run.  Comparisons  of  tide  and  velocity  data  of  RMA-2V  inset 
mesh,  RMA-2V  global  mesh,  and  physical  model  were  good.  The  comparisons  of 
tide  and  velocity  data  for  the  inset  and  global  mesh  stations  showed  no  gen¬ 
eral  pattern  of  variation,  phase  agreement  was  good,  and  tide  and  velocity 
amplitudes  differed  from  0  to  0.5  ft  (tide)  and  0  to  0.5  fps  (velocity).  When 
the  inset  mesh  and  physical  model  data  were  compared,  the  RMA-2V  velocity  val¬ 
ues  again  tended  to  be  slightly  lower  than  the  physical  model  values.  How¬ 
ever,  the  inset  mesh  velocity  comparisons  were  slightly  better  than  global 
mesh  versus  physical  model  comparisons.  This  was  due  to  the  higher  resolution 
of  the  smaller  inset  mesh.  Plate  15  is  the  plotted  data  from  selected  sta¬ 
tions  of  inset  mesh  versus  global  mesh,  and  Plate  16  shows  data  for  the  inset 
mesh  versus  the  physical  model. 

Results 

A5.  The  RMA-2V  hydrodynamic  model  was  verified  to  the  two  sets  of  the 
physical  model  data.  Both  the  global  and  the  high-resolution  inset  mesh  re¬ 
produced  the  tides  and  velocities  of  the  SFBM  with  reasonable  accuracy  and 
repeatability . 


Disposal  Model  Validation 


A6.  In  this  application  the  disposal  model,  DIFID,  was  not  verified  to 
field  or  other  model  data.  It  has,  however,  been  shown  to  provide  reasonable 
results  at  this  site  compared  to  field  experience  and  good  comparison  to  field 
data  at  other  sites  (Trawle  and  Johnson  1986b).  Thus  the  disposal  model  is 
considered  to  be  validated  (well  grounded  on  principles  and  generally  capable 
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of  reproducing  observed  behavior)  but  not  verified  (compared  directly  to  San 
Francisco  Bay/Alcatraz  site  field  data). 

Test  conditions 

47.  The  disposal  simulation  consisted  of  1,000  cu  yd  of  hydraulically 
dredged  material,  primarily  silt-clays,  being  disposed  from  a  stationary  barge 
in  the  Alcatraz  dredged  material  disposal  site.  Sediment  characteristics  and 
barge  location  were  similar  to  the  information  presented  in  Phase  II  of  the 
Alcatraz  disposal  site  survey  (Science  Applications  International  Corporation 
(SAIC)  1987). 

48.  Two  disposal  operations  were  simulated,  one  at  maximum  ebb  and  the 
other  at  maximum  flood  velocities.  The  clay-silt  material  that  did  not  reach 
bottom  but  remained  in  suspension  was  tabulated  at  four  depths,  10,  30,  50, 
and  70  ft,  at  several  time  intervals  after  the  disposal  operation. 

49.  The  DIFID  model  was  run  using  a  velocity  field  from  the  RMA-2V  mean 
tide  simulation  (Figure  9).  A  special  version  of  VPLOT  converted  the  velocity 
vectors  from  the  nodes  and  midsides  of  the  elements  in  the  inset  mesh  to  a 
regular  mesh  of  the  size  and  spacing  designed  (Figure  10). 

Model  coefficients 

50.  There  are  14  coefficients  in  DIFID  that  require  input  values.  De¬ 
fault  values  based  on  previous  laboratory  and  field  experiments  are  included 
in  the  model.  Computer  experimentation  such  as  presented  by  Johnson  and 
Holliday  (1978)  has  shown  that  model  results  appear  to  be  fairly  insensitive 
to  many  of  the  coefficients.  Therefore,  the  model  default  values  were  used 
for  these  simulations. 

Disposal  model  results 

51.  The  suspended  sediment  plume  predicted  by  DIFID  displayed  similar 
characteristics  to  the  SAIC  (1987)  Phase  II  study  in  that  it  did  not  persist 
beyond  10-15  min,  had  comparable  flood  and  ebb  plume  direction  and  distance  of 
travel  features,  and  had  no  significant  north/south  dispersion.  Figure  11  is 
a  sample  of  the  printed  output  from  the  DIFID  ebb  velocity  simulation  showing 
silt-clay  concentrations  40  ft  below  the  water  surface,  540  sec  (9  min)  after 
the  disposal.  Alcatraz  Island  (land)  and  the  location  of  the  disposal  barge 
are  shown  on  the  grid. 

Discussion  of  results 

52.  The  field  data  available  in  the  SAIC  (1987)  report  indicate  the 
disposal  model  reasonably  simulated  the  disposal  process  in  the  Alcatraz  site. 
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Alcatraz  Island,  RMA-2V  velocity  vectors  at  maximum  ebb 
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Sample  DIFID  output,  summary  of:  silt-clay  distribution  after  540  sec,  40  ft 

below  water  surface 


The  default  model  coefficients  used  in  DIFID  could  be  made  more  site  specific 
in  the  future  as  a  better  understanding  of  sediment  behavior  at  the  Alcatraz 
site  is  achieved. 

53.  The  purpose  of  the  simulation  was  to  produce  disposal-related  sus¬ 
pended  sediment  concentrations  as  input  into  the  TABS-2  sediment  transport 
model  STUDH.  This  was  achieved  by  averaging  the  suspended  sediment  concentra¬ 
tions  in  the  water  column  in  each  cell  of  the  DIFID  orid  and  assigning  these 
values  to  the  comparable  inset  mesh  nodes  and  midsides. 

Application  of  Sediment  Transport  Model 


54.  The  sediment  transport  model  STUDH  was  used  to  trace  the  sediment 
movement  through  the  system.  As  in  the  case  of  the  disposal  model,  no  field 
data  were  available  to  verify  STUDH.  Normally  in  the  verification  process, 
the  model  coefficients  are  adjusted  until  the  model  produces  results  similar 
to  conditions  observed  in  the  prototype.  In  the  absence  of  verification  data, 
the  coefficients  used  were  realistic  estimates  of  the  general  prototype  and 
sediment  characteristics  found  in  the  Alcatraz  disposal  area. 

Test  conditions 

55.  The  inset  mesh  geometry  was  used  with  the  RMA-2V  hydrodynamic  out¬ 
put  data  from  the  mean  tide,  low  flow  simulation.  Along  with  the  run  control 
file,  these  data,  converted  into  metric  units,  are  necessary  input  to  run 
STUDH.  Sediment  concentrations  (from  DIFID)  assigned  to  designated  inset  mesh 
nodes  and  the  sediment  transport  coefficients  were  included  in  the  run  control 
file.  The  sediment  concentrations  from  the  disposal  at  maximum  ebb  velocity 
were  used,  and  seven  time-steps  (two  before  and  four  after  maximum  ebb  veloc¬ 
ity  for  a  total  of  4.3  hr)  were  simulated. 

Model  coefficients 

56.  The  primary  function  of  this  simulation  was  to  evaluate  the  use  of 
STUDH  in  the  application  of  tracing  suspended  sediment  movement  in  the  system. 
Therefore,  all  terms  relating  to  bed  deposition  and  erosion  were  turned  off. 
Only  required  coefficients  were  used.  A  fluid  temperature  of  20°  C  and  a 
fluid  density  of  1,025  kg/cu  m  were  also  specified. 

57.  The  initial  runs  of  STUDH  were  unstable  with  oscillations  in  the 
sediment  concentration  results.  Adjustments  were  made  one  at  a  time  to  the 
length  '•f  th*1  tim*»-etep  (37  min,  18  min,  and  4  min);  x-  and  y-diffusion 
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coefficients  (EX  =  25  sq  m/sec,  EY  =  10  sq  m/sec;  and  EX  =  50  sq  m/sec,  EY  = 

25  sq  m/sec);  and  Crank-Nicholson  THETA  values  (0.5,  0.67,  0.9,  and  1.0). 
Simulations  were  also  made  with  and  without  background  and  boundary  sediment 
concentrations.  The  adjustment  that  made  the  most  difference  and  improved  the 
stability  of  the  model  was  the  Crank-Nicholson  THETA  value.  With  a  value  of 
1.0,  the  oscillations  were  quelled  and  the  solutions  were  stable.  The  simula¬ 
tion  described  in  this  report  was  run  with  the  following:  time-step  of 
2,235.6  sec  (37  min);  diffusion  coefficients  of  EX  25  sq  m/s ec ,  EY  10  sq 
m/sec;  Crank-Nicholson  THETA  of  1.0;  and  a  background  and  boundary  sediment 
concentration  of  20  mg/i. 

Transport  model  results 

58.  The  STUDH  results,  although  unverified,  were  promising  and  showed  a 
reasonable  progression  of  suspended  sediment  out  of  the  disposal  site.  The 
results  were  more  stable  when  the  simulation  was  run  with  a  background  concen¬ 
tration,  the  solution  ending  at  the  background  value.  When  no  background  was 
specified,  the  solution  continued  to  approach  zero  but  never  quite  reached  it, 
sometimes  going  to  small  negative  numbers.  Plates  17-19  are  plots  of  the 
predicted  suspended  sediment  at  37-min  time-steps  for  3.73  hr  after  the  dis¬ 
posal  operation,  which  took  place  during  the  time  of  ebb  current  velocity. 
Results 

59.  It  has  been  demonstrated  that  STUDH  can  be  used  to  trace  the  move¬ 
ment  of  suspended  sediment  from  the  Alcatraz  disposal  site  through  the  Central 
Bay  system.  Additional  refinement  and  prototype  data  are  needed  to  adjust  the 
model  coefficients  to  be  more  representative  of  the  transport  process. 

60.  In  the  absence  of  verification  field  data,  the  modelers  must  rely 
on  their  experience  and  familiarity  with  both  the  San  Francisco  Bay  and  the 
TABS-2  modeling  systems  to  establish  the  modeling  coefficients  that  best 
represent  the  prototype  processes. 
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PART  IV:  CONCLUSIONS 


61.  The  use  of  the  San  Francisco  Bay-Delta  (SFBM)  physical  model  and 
the  TABS-2  modeling  system  has  produced  an  interrelated  set  of  models  that  can 
simulate  and  address  many  hydrodynamic,  circulation,  and  suspended  sediment 
transport  questions  of  most  portions  of  the  San  Francisco  Bay  system. 

62.  The  RMA-2V  hydrodynamic  model  of  both  the  comprehensive  and  the 
more  detailed  inset  meshes  have  been  verified  to  data  produced  in  the  SFBM 
physical  model. 

63.  The  disposal  model,  DIFID,  produced  suspended  sediment  results 
similar  to  field  data  and  established  a  realistic  suspended  sediment  cloud 
that  might  be  produced  by  a  dredged  material  disposal  operation  in  the 
Alcatraz  disposal  site  during  the  time  following  maximum  ebb  current 
velocities . 

64.  The  coupling  of  the  geometry,  hydrodynamics,  and  sediment  input  ’n 
the  transport  model,  STUDH,  produced  a  reasonable  sediment  transport  pattern 
out  of  the  disposal  site.  Specific  clouds  of  sediment  concentrations  could  be 
identified  and  monitored  as  they  travelled  and  dispersed  through  the  system. 

65.  The  modeling  system  has  its  capabilities  and  applications.  How¬ 
ever,  the  results  are  just  reasonable  simulations,  not  fully  verified  ones. 
Each  of  the  numerical  models,  RMA-2V,  DIFID,  and  STUDH,  has  individual  capa¬ 
bilities  and  limitations,  the  greatest  of  which  is  the  two-dimensional  approx¬ 
imation  of  a  three-dimensional  phenomenon.  The  vertically  averaged  velocities 
and  sediment  fields  will  mask  two-layer  flow  and  other  three-dimensional  pro¬ 
cesses.  Even  with  this  simplification,  the  model  results  are  useful  in  esti¬ 
mating  the  short-  and  long-term  fates  of  sediments  released  during  a  disposal 
operation. 
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APPENDIX  A:  THE  TABS-2  SYSTEM 


1.  TABS-2  is  a  collection  of  generalized  computer  programs  and  utility 
codes  integrated  into  a  numerical  modeling  system  for  studying  two-dimensional 
hydrodynamics,  sedimentation,  and  transport  problems  in  rivers,  reservoirs, 
bays,  and  estuaries.  A  schematic  representation  of  the  system  is  shown  in 
Figure  Al.  It  can  be  used  either  as  a  stand-alone  solution  technique  or  as  a 
step  in  the  hybrid  modeling  approach.  The  basic  concept  is  to  calculate 
water-surface  elevations,  current  patterns,  sediment  erosion,  transport  and 
deposition,  the  resulting  bed  surface  elevations,  and  the  feedback  to  hydrau¬ 
lics.  Existing  and  proposed  geometry  can  be  analyzed  to  determine  the  impact 
on  sedimentation  of  project  designs  and  to  determine  the  impact  of  project 
designs  on  salinity  and  on  the  stream  system.  The  system  is  described  in  de¬ 
tail  by  Thomas  and  McAnally  (1985). 

2.  The  three  basic  components  of  the  system  are  as  follows: 

a.  "A  Two-Dimensional  Model  for  Free  Surface  Flows,"  RMA-2V, 

b.  "Sediment  Transport  in  Unsteady  2-Dimensional  Flows,  Horizontal 
Plane,"  STUDH. 

c.  "Two-Dimensional  Finite  Element  Program  for  Water  Quality," 

RMA-4 . 

3.  RMA-2V  is  a  finite  element  solution  of  the  Reynolds  form  of  the 
Navier-Stokes  equations  for  turbulent  flows.  Friction  is  calculated  with 
Manning's  equation  and  eddy  viscosity  coefficients  are  used  to  define  the 
turbulent  losses.  A  velocity  form  of  the  basic  equation  is  used  with  side 
boundaries  treated  as  either  slip  or  static.  The  model  automatical lv  recog¬ 
nizes  dry  elements  and  corrects  the  mesh  accordingly.  Boundary  conditions  may 
be  water-surface  elevations,  velocities,  or  discharges  and  may  occur  inside 
the  mesh  as  well  as  along  the  edges. 


TABS-2 


Figure  Al.  TABS-2  schematic 
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4.  The  sedimentation  model,  STUDH,  solves  the  convection-diffusion 
equation  with  bed  source  terms.  These  terms  are  structured  for  either  sand  or 
cohesive  sediments.  The  Ackers-White  (1973)  procedure  is  used  to  calculate  a 
sediment  transport  potential  for  the  sands  from  which  the  actual  transport  is 
calculated  based  on  availability.  Clay  erosion  is  based  on  work  by  Parthen- 
iades  (1962)  and  Ariathurai  and  the  deposition  of  clay  utilizes  Krone's  equa¬ 
tions  (Ariathurai,  MacArthur,  and  Krone  1977).  Deposited  material  forms 
layers,  as  shown  in  Figure  A2,  and  bookkeeping  allows  up  to  10  layers  at  each 
node  for  maintaining  separate  material  types,  deposit  thickness,  and  age.  The 
code  uses  the  same  mesh  as  RMA-2V. 

5.  Salinity  calculations,  RMA-4 ,  are  made  with  a  form  of  the 
convective-diffusion  equation  which  has  general  source-sink  terms.  Up  to 
seven  conservative  substances  or  substances  requiring  a  decay  term  can  be 
routed.  The  code  uses  the  same  mesh  as  RMA-2V. 

6.  Each  of  these  generalized  computer  codes  can  be  used  as  a  stand¬ 
alone  program,  but  to  facilitate  the  preparation  of  input  data  and  to  aid  in 
analyzing  results,  a  family  of  utility  programs  was  developed  for  the  follow¬ 
ing  purposes: 


a . 

Digitizing 

b. 

Mesh  generation 

c . 

Spatial  data  management 

d. 

Graphical  output 

e . 

Output  analysis 

f . 

File  management 

g- 

Interfaces 

h. 

Job  control  language 

Finite  Element  Modeling 

7.  The  TABS-2  numerical  models  used  in  this  effort  employ  the  finite 
element  method  to  solve  the  governing  equations.  To  help  those  who  are  un¬ 
familiar  with  the  method  to  better  understand  this  report,  a  brief  description 
of  the  method  is  given  here. 

8.  The  finite  element  method  approximates  a  solution  to  equations  by 
dividing  the  area  of  interest  into  smaller  subareas,  which  are  called  ele¬ 
ments.  The  dependent  variables  (e.g.,  water-surface  elevations  and  sediment 
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concentrations)  are  approximated  over  each  element  by  continuous  functions 
which  interpolate  in  terms  of  unknown  point  (node)  values  of  the  variables. 

An  error,  defined  as  the  deviation  of  the  approximation  solution  from  the  cor¬ 
rect  solution,  is  minimized.  Then,  when  boundary  conditions  are  imposed,  a 
set  of  solvable  simultaneous  equations  is  created.  The  solution  is  continuous 
over  the  area  of  interest. 

9.  In  one-dimensional  problems,  elements  are  line  segments.  In  two- 
dimensional  problems,  the  elements  are  polygons,  usually  either  triangles  or 
quadrilaterals.  Nodes  are  located  on  the  edges  of  elements  and  occasionally 
inside  the  elements.  The  interpolating  functions  may  be  linear  or  higher 
order  polynomials.  Figure  A2  illustrates  a  quadrilateral  element  with  eight 
nodes  and  a  linear  solution  surface  where  F  is  the  interpolating  function. 

10.  Most  water  resource  applications  of  the  finite  element  method  use 
the  Galerkin  method  of  weighted  residuals  to  minimize  error.  In  this  method 
the  residual,  the  total  error  between  the  approximate  and  correct  solutions, 
is  weighted  by  a  function  that  is  identical  with  the  interpolating  function 
and  then  minimized.  Minimization  results  in  a  set  of  simultaneous  equations 
in  terms  of  nodal  values  of  the  dependent  variable  (e.g.  water-surface  eleva¬ 
tions  or  sediment  concentration).  The  time  portion  of  time-dependent  problems 
can  be  solved  by  the  finite  element  method,  but  it  is  generally  more  efficient 
to  express  derivatives  with  respect  to  time  in  finite  difference  form. 

The  Hydrodynamic  Model,  RMA-2V 


Applications 

11.  This  program  is  designed  for  far-field  problems  in  which  vertical 
accelerations  are  negligible  and  the  velocity  vectors  at  a  node  generally 
point  in  the  same  directions  over  the  entire  depth  of  the  water  column  at  any 
Instant  of  time.  It  expects  a  homogeneous  fluid  with  a  free  surface.  Both 
steady  and  unsteady  state  problems  can  be  analyzed.  A  surface  wind  stress  can 
be  imposed. 

12.  The  program  has  been  applied  to  calculate  flow  distribution  around 
islands;  flow  at  bridges  having  one  or  more  relief  openings,  in  contracting 
and  expanding  reaches,  into  and  out  of  off-channel  hydropower  plants,  at  river 
junctions,  and  into  and  out  of  pumping  plant  channels;  and  general  flow  pat¬ 
terns  in  rivers,  reservoirs,  and  estuaries. 
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Limitations 


13.  This  program  is  not  designed  for  near-field  problems  where  flow- 
structure  interactions  (such  as  vortices,  vibrations,  or  vertical  accelera¬ 
tions)  are  of  interest.  Areas  of  vertically  stratified  flow  are  beyond  this 
program's  capability  unless  it  is  used  in  a  hybrid  modeling  approach.  It  is 
two-dimensional  in  the  horizontal  plane,  and  zones  where  the  bottom  current  is 
in  a  different  direction  from  the  surface  current  must  be  analyzed  with  con¬ 
siderable  subjective  judgement  regarding  long-term  energy  considerations.  It 
is  a  free-surface  calculation  for  subcritical  flow  problems. 

Governing  equations 

14.  The  generalized  computer  program  RMA-2V  solves  the  depth- integrated 
equations  of  fluid  mass  and  momentum  conservation  in  two  horizontal  direc¬ 
tions.  The  form  of  the  solved  equations  is 
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e  =  eddy  viscosity  coefficient,  for  xx  =  normal  direction  on 

x-axis  surface;  yy  =  normal  direction  on  y-axis  surface;  xy 
and  yx  =  shear  direction  on  each  surface 

g  =  acceleration  due  to  gravity 

a  =  elevation  of  bottom 

n  =  Manning's  n  value 

1.486  =  conversion  from  SI  (metric)  to  non-SI  units 
5  =  empirical  wind  shear  coefficient 

V  =  wind  speed 
a 

ip  =  wind  direction 

w  =  rate  of  earth's  angular  rotation 
cf  =  local  latitude 

15.  Equations  A1 ,  A2,  and  A3  are  solved  by  the  finite  element  method 
using  Galcvkin  weighted  residuals.  The  elements  may  be  either  quadrilaterals 
or  triangles  and  may  have  curved  (parabolic)  sides.  The  shape  functions  are 
quadratic  for  flow  and  linear  for  depth.  Integration  in  space  is  performed  by 
Gaussian  integration.  Derivatives  in  time  are  replaced  by  a  nonlinear  finite 
difference  approximation.  Variables  are  assumed  to  vary  over  each  time  inter¬ 
val  in  the  form 

f(t)  =  f(0)  +  at  +  btC  tQ  <  t  <  t  (A4) 

which  is  differentiated  with  respect  to  time,  and  cast  in  finite  difference 
form.  Letters  a  ,  b  ,  and  c  are  constants.  It  has  been  found  by  experi¬ 
ment  that  the  best  value  for  c  is  1.5  (Norton  and  King  1977). 

16.  The  solution  is  fully  implicit  and  the  set  of  simultaneous  equations 
is  solved  by  Newton-Raphson  iteration.  The  computer  code  executes  the  solu¬ 
tion  by  means  of  a  front-type  solver  that  assembles  a  portion  of  the  matrix 
and  solves  it  before  assembling  the  next  portion  of  the  matrix.  The  front 
solver's  efficiency  is  largely  independent  of  bandwidth  and  thus  does  not 
require  as  much  care  in  formation  of  the  computational  mesh  as  do  traditional 
solvers . 

17.  The  code  RMA-2V  is  based  on  the  earlier  version  RMA-2  (Norton  and 
King  1977)  but  differs  from  it  in  several  ways.  It  is  formulated  in  terms  of 
velocity  (v)  instead  of  unit  discharge  (vh) ,  which  improves  some  aspects  of 
the  code's  behavior;  it  permits  drying  and  wetting  of  areas  within  the  grid; 
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and  it  permits  specification  of  turbulent  exchange  coefficients  in  directions 
other  than  along  the  x-  and  z-axes.  For  a  more  complete  description,  see 
Appendix  F  of  Thomas  and  McAnally  (1985). 

The  Sediment  Transport  Model,  STUDH 


Applications 

18.  STUDH  can  be  applied  to  clay  and/or  sand  bed  sediments  where  flow 
velocities  can  be  considered  two-dimensional  (i.e.,  the  speed  and  direction 
can  be  satisfactorily  represented  as  a  depth-averaged  velocity).  If  is  useful 
for  both  deposition  and  erosion  studies  and,  to  a  limited  extent,  for  stream 
width  studies.  The  program  treats  two  categories  of  sediment:  noncohesive, 
which  is  referred  to  as  sand  here,  and  cohesive,  which  is  referred  to  as  clay. 
Limitations 

19.  Both  clay  and  sand  may  be  analyzed,  but  the  model  considers  a 
single,  effective  grain  size  for  each  and  treats  each  separately.  Fall  veloc¬ 
ity  must  be  prescribed  along  with  the  water-surface  elevations,  x-velocity, 
y-velocity,  diffusion  coefficients,  bed  density,  critical  shear  stresses  for 
erosion,  erosion  ra'_e  constants,  and  critical  shear  stress  for  deposition. 

20.  Many  applications  cannot  use  long  simulation  periods  because  of 
their  computation  cost.  Study  areas  should  be  made  as  small  as  possible  to 
avoid  an  excessive  number  of  elements  when  dynamic  runs  are  contemplated  yet 
must  be  large  enough  to  permit  proper  posing  of  boundary  conditions.  The  same 
computation  time  interval  must  be  satisfactory  for  both  the  transverse  and 
longitudinal  flow  directions. 

21.  The  program  does  not  compute  water-surface  elevations  or  velocities; 
therefore  these  data  must  be  provided.  For  complicated  geometries,  the  numer¬ 
ical  model  for  hydrodynamic  computations,  RMA-2V,  is  used. 

Governing  equations 

22.  The  generalized  computer  program  STUDH  solves  the  depth-integrated 
convection-dispersion  equation  in  two  horizontal  dimensions  for  a  single  sedi¬ 
ment  constituent.  For  a  more  complete  description,  see  Appendix  G  of  Thomas 
and  McAnally  (1985).  The  form  cf  the  solved  equation  is 
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where 


C  =  concentration  of  sediment 
u  =  depth-integrated  velocity  in  x-direction 
v  =  depth-integrated  velocity  in  y-direction 
D  =  dispersion  coefficient  in  x-direction 
D  =  dispersion  coefficient  in  y-direction 

=  coefficient  of  concentration-dependent  source/sink  term 
a ^  =  coefficient  of  source/sink  term 

23.  The  source/sink  terms  in  Equation  A5  are  computed  in  routines  that 
treat  the  interaction  of  the  flow  and  the  bed.  Separate  sections  of  the  code 
handle  computations  for  clay  bed  and  sand  bed  problems. 

Sand  transport 

24.  The  source/sink  terms  are  evaluated  by  first  computing  a  potential 
sand  transport  capacity  for  the  specified  flow  conditions,  comparing  that 
capacity  with  the  amount  of  sand  actually  being  transported,  and  then  eroding 
from  or  depositing  to  the  bed  at  a  rate  that  would  approach  the  equilibrium 
value  after  sufficient  elapsed  time. 

25.  The  potential  sand  transport  capacity  in  the  model  is  computed  by 
the  method  of  Ackers  and  White  (1973)  ,  which  uses  a  transport  power  (work 
rate)  approach.  It  has  been  shown  to  provide  superior  results  for  transport 
under  steady-flow  conditions  (White,  Milli,  and  Crabbe  1975)  and  for  combined 
waves  and  currents  (Swart  1976).  Flume  tests  at  the  US  Army  Engineer  Water¬ 
ways  Experiment  Station  have  shown  that  the  concept  is  valid  for  transport  by 
estuarine  currents. 

26.  The  total  load  transport  function  of  Ackers  and  White  is  based  upon 
a  dimensionless  grain  size 
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where 

D  =  sediment  particle  diameter 
s  =  specific  gravity  of  the  sediment 
v  =  kinematic  viscosity  of  the  fluid 
and  a  sediment  mobility  parameter 
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F 

gr 

where 

x  =  total  boundary  shear  stress 

n'  =  a  coefficient  expressing  the  relative  importance  of  bed-load  and 
suspended-load  transport,  given  in  Equation  A9 

t'  =  boundary  surface  shear  stress 

The  surface  shear  stress  is  that  part  of  the  total  shear  stress  which  is  due 
to  the  rough  surface  of  the  bed  only,  i.e.,  not  including  that  part  due  to  bed 
forms  and  geometry.  It  therefore  corresponds  to  that  shear  stress  that  the 
flow  would  exert  on  a  plane  bed. 

27.  The  total  sediment  transport  is  expressed  as  an  effective 
concentration 
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For  D  <60 
gr 


n*  =  0.00 


( A 1  3  ) 
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A  =  0.17 


(A14) 


C  =  0.025  (A15) 

m  =  1.5  (A16) 

28.  Equations  A6-A16  result  in  a  potential  sediment  concentration  . 
This  value  is  the  depth-averaged  concentration  of  sediment  that  will  occur  if 
an  equilibrium  transport  rate  is  reached  with  a  nonlimited  supply  of  sediment. 
The  rate  of  sediment  deposition  (or  erosion)  is  then  computed  as 


where 


C  =  present  sediment  concentration 

t  =  time  constant 
c 

For  deposition,  the  time  constanc  is 
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and  for  erosion  it  is 
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where 

At  =  computational  time-step 

=  response  time  coefficient  for  deposition 

=  sediment  settling  velocity 

C  =  response  time  coefficient  for  erosion 
e 
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The  sand  bed  has  a  specified  initial  thickness  which  limits  the  amount  of  ero¬ 
sion  to  that  thickness. 

Cohesive  sediments  transport 

29.  Cohesive  sediments  (usually  clays  and  some  silts)  are  considered  to 
be  depositional  if  the  bed  shear  stress  exerted  by  the  flow  is  less  than  a 
critical  value  .  When  that  value  occurs,  the  deposition  rate  is  given  by 

Krone's  (1962)  equation 
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where 

S  =  source  term 

=  fall  velocity  of  a  sediment  particle 
h  =  flow  depth 

C  =  sediment  concentration  in  water  column 
t  =  bed  shear  stress 

=  critical  shear  stress  for  deposition 
C^  =  critical  concentration  =  300  mg/i 


30.  If  the  bed  shear  stress  is  greater  than  the  critical  value  for  par¬ 
ticle  erosion  t  ,  material  is  removed  from  the  bed.  The  source  term  is  then 
e 

computed  by  Ariathurai's  (Ariathurai,  MacArthur,  and  Krone  1977)  adaptation  of 
Partheniades '  (1962)  findings: 
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where  P  is  the  erosion  rate  constant,  unless  the  shear  stress  is  also 
greater  than  the  critical  value  for  mass  erosion.  When  this  value  is 
exceeded,  mass  failure  of  a  sediment  layer  occurs  and 
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where 

=  th'ckness  of  the  failed  layer 

=  density  of  the  failed  layer 

At  =  time  interval  over  which  failure  occurs 

r  =  bulk  shear  strength  of  the  layer 
s 

31,  The  cohesive  sediment  bed  consists  of  1  to  10  layers,  each  with  a 
distinct  density  and  erosion  resistance.  The  layers  consolidate  with 
overburden  and  time. 

Bed  shear  stress 

32.  Bed  shear  stresses  are  calculated  from  the  flow  speed  according  to 
one  of  four  optional  equations:  the  smooth-wall  log  velocity  profile  or 
Manning  equation  for  flows  alone;  and  a  smooth  bed  ut  rippled  bed  equation  for 
combined  currents  and  wind  waves.  Shear  stresses  are  calculated  using  the 
shear  velocity  concept  where 
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where 

t,  =  bed  shear  stress 
b 

u^  =  shear  velocity 

and  the  shear  velocity  is  calculated  by  one  of  four  methods: 
a.  Smooth-wall  log  velocity  profiles 
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which  is  applicable  to  the  lower  15  percent  of  the  boundary 
layer  when 
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t>  30 
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where  u  is  the  mean  flow  velocity  (resultant  of  u  and  v 
components) 

b.  The  Manning  shear  stress  equation 
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where  CME  is  a  coefficient  of  1  for  SI  (metric)  units  and 
1.486  for  non-SI  units  of  measurement. 

c.  A  Jonsson-type  equation  for'  surface  shear  stress  (plane  beds) 
caused  by  waves  and  currents 
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where 


f  =  shear  stress  coefficient  for  waves 
w 

u  =  maximum  orbital  velocity  of  waves 
om 

f  =  shear  stress  coefficient  for  currents 
c 


d,  A  Bijker-type  equation  for  total  shear  ..tress  caused  by  waves 
and  current 
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solution  method 

13.  Equation  A5  is  solved  by  the  finite  element  method  using  Caierkin 
weighted  residuals.  Like  RMA-2V,  which  uses  the  same  general  solution  tech¬ 
nique,  elements  are  quadrilateral  and  may  have  parabolic  sides.  Shape  func¬ 
tions  are  quadratic.  Integration  in  space  is  Gaussian.  Time-stepping  is 
performed  by  a  Crank-Ni cholson  approach  with  a  weighting  factor  (0)  of  0 .bo. 

A13 


A  front-type  solver  similar  to  that  in  RMA-2V  is  used  to  solve  the  simultane¬ 
ous  equations. 
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APPENDIX  B:  DESCRIPTION  OF  THE  NUMERICAL  MODEL,  DIFID* 


Model  Origin 

1.  The  instantaneous  dump  model  (DIFID)  was  developed  by  Brandsma  and 
Divoky  (1976)  for  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  un¬ 
der  the  Dredged  Material  Research  Program.  Much  of  the  basis  for  the  model 
was  provided  by  earlier  model  development  by  Koh  and  Chang  (1973)  for  the 
barged  disposal  of  wastes  in  the  ocean.  That  work  was  conducted  under  funding 
by  the  Environmental  Protection  Agency  in  Corvallis,  Oregon.  Modifications  to 
the  original  model  have  been  made  by  the  Hydraulics  Laboratory  at  WES. 

Model  Approach 


2.  The  model  simulates  movement  of  the  disposed  material  as  it  falls 
through  the  water  column,  spreads  over  the  bottom,  and  finally  is  transported 
and  diffused  as  suspended  sediment  by  the  ambient  current.  DIFID  is  designed 
to  simulate  the  movement  of  material  from  an  instantaneous  dump  which  falls  at 
a  hemispherical  cloud.  Thus  the  total  time  required  for  the  material  to  leave 
the  disposal  vessel  should  not  be  substantially  greater  than  the  time  required 
for  the  material  to  reach  the  bottom. 

3.  The  model  requires  that  the  dredged  material  be  broken  into  various 
solid  fractions  with  a  settling  velocity  specified  for  each  fraction.  In  many 
cases,  a  significant  portion  of  the  material  falls  as  "clumps"  which  may  have 
a  settling  velocity  of  perhaps  1.0  to  5.0  fps.  This  is  especially  true  if  the 
dredging  is  done  by  clamshell  and  can  be  true  in  the  case  ef  hydraulicall v 
dredged  material  if  consolidation  takes  place  in  the  hopper  during  transit  to 
the  disposal  site.  The  specification  of  a  "clump"  fraction  is  rather  subjec¬ 
tive;  therefore  the  inability  to  accurately  characterize  the  disposed  material 
in  some  disposal  operations  prevents  a  quantitative  interpretation  of  model 
results  in  those  operations. 

4.  As  noted,  a  settling  velocity  must  be  prescribed  for  each  solid 
fraction.  A  basic  assumption  is  that  unless  the  fraction  is  specified  as  be¬ 
ing  cohesive,  in  which  case  the  settling  velocity  is  computed  as  a  function  of 


*  Taken  from  Trawle  and  Johnson  (1986),  pp  7-18. 
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concentration,  the  settling  is  considered  to  occur  at  a  constant  rate.  In 
other  words,  hindered  settling  is  not  taken  into  consideration. 

5.  Although  a  variable  water  depth  is  allowed  over  the  long-term  grid, 
the  collapse  of  the  dredged  material  cloud  on  the  bottom  is  somewhat  re¬ 
stricted.  The  effect  of  a  bottom  slope  is  allowed  through  the  incorporation 
of  a  gravitational  force  in  the  computation  of  the  collapsing  cloud.  However, 
a  basic  limitation  still  exists  in  that  the  bottom  is  assumed  to  slope  in  only- 
one  direction  over  the  collapsed  region;  e.g.,  bottom  collapse  on  a  "mound” 
where  the  collapsing  cloud  runs  down  the  sides  is  not  treated. 

6.  A  major  limitation  of  the  model  is  the  basic  assumption  that  once 
solid  particles  are  deposited  on  the  bottom  they  remain  there.  Therefore  the 
model  should  only  be  applied  over  time  frames  in  which  erosion  of  the  newly 
deposited  material  is  insignificant. 

7.  The  model  allows  for  two  separate  treatments  of  the  passive  trans¬ 
port  and  diffusion  phase.  In  one  method,  material  from  the  convective  descent 
and  dynamic  collapse  phases  is  inserted  into  the  fixed  long-term  grid.  There¬ 
fore  computations  at  each  point  of  the  grid  must  be  made  at  each  time-step  in 
order  to  march  the  solution  from  one  time-step  to  the  next.  Solid  bodies  in 
the  field  and  boundary  effects  are  treated.  However,  disadvantages  are  that 
the  vertical  distribution  is  assumed  to  be  that  of  a  "top  hat"  profile,  grid 
dispersion  errors  may  occur,  and  the  computations  can  become  costly  for  large 
grids  if  many  time-steps  are  computed.  A  "top  hat"  profile  is  represented  by 
a  step  function  that  does  not  allow  for  any  gradual  change  over  the  water  col¬ 
umn.  The  second  approach  is  to  allow  material  from  the  descent  and  collapse 
phases  to  be  stored  in  small  Gaussian  clouds.  These  clouds  are  then  diffused 
and  transported  at  the  end  of  each  time-step.  Computations  on  the  long-term 
grid  are  only  made  at  those  times  when  output  is  desired.  However,  a  limita¬ 
tion  when  using  this  approach  is  that  horizontal  solid  boundaries  are  not 
allowed  in  the  long-term  grid.  This  limitation  could  be  removed  by  employing 
reflection  principles  as  is  currently  done  at  the  surface  and  the  bottom. 

Theoretical  Basis 


8.  The  behavior  of  the  disposed  material  is  assumed  to  be  separated 
into  three  phases:  convective  descent,  during  which  the  dump  cloud  or  dis¬ 
charge  jet  falls  under  the  influence  of  gravity;  dynamic  collapse,  occurring 
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when  the  descending  cloud  either  impacts  the  bottom  or  arrives  at  the  level  of 
neutral  buoyancy  at  which  descent  is  retarded  and  horizontal  spreading  domi¬ 
nates;  and  long-term  passive  dispersion,  commencing  when  the  material  trans¬ 
port  and  spreading  are  determined  more  by  ambient  currents  and  turbulence  than 
by  the  dynamics  of  the  disposal  operation.  Figure  Bl  illustrates  these 
phases . 


Figure  Bl.  Illustration  of  idealized  bottom  encounter  after  instantaneous 
dump  of  dredged  material  (from  Brandsma  and  Divoky  1976) 


Convective  descent 

9.  A  single  cloud  that  maintains  a  hemispherical  shape  during  convec¬ 
tive  descent  is  assumed  to  be  released.  Since  the  solids  concentration  in 
discharged  dredged  material  is  usually  low,  the  cloud  is  expected  to  behave  as 
a  dense  liquid;  thus  a  basic  assumption  is  that  a  buoyant  thermal  analysis  is 
appropriate.  The  equations  governing  the  motion  are  those  for  conservation  of 
mass,  momentum,  buoyancy,  each  solid,  and  vorticity.  The  equations  are 
straightforward  statements  of  conservation  principles  and  are  presented  by 
Brandsma  and  Divoky  (1976).  It  should  be  noted  that  the  entrainment  coeffi¬ 
cient  associated  with  the  entrainment  of  ambient  fluid  into  the  descending 
hemispherical  cloud  is  assumed  to  vary  smoothly  between  its  value  for  a  vertex 


ring  and  the  value  for  turbulent  thermals.  Model  results  are  quite  sensitive 
to  the  entrainment  coefficient,  which  in  turn  is  dependent  upon  the  material 
being  dumped  (the  higher  the  moisture  content,  the  larger  the  value  of  the 
entrainment  coefficient). 

Dynamic  collapse 

10.  During  convective  descent,  the  dumped  material  cloud  grows  as  a  re¬ 
sult  of  entrainment.  Eventually,  either  the  material  reaches  the  bottom  or 
the  density  difference  between  the  discharged  material  and  the  ambient  fluid 
becomes  small  enough  for  a  position  of  neutral  buoyancy  to  be  assumed.  In 
either  case,  the  vertical  motion  is  arrested  and  a  dynamic  spreading  in  the 
horizontal  direction  occurs.  The  basic  shape  assumed  for  the  collapsing  cloud 
is  an  oblate  spheroid.  With  the  exception  of  vorticity,  which  is  assumed  to 
have  been  dissipated  by  the  stratified  ambient  water  column,  the  same  conser¬ 
vation  equations  used  in  convective  descent  but  now  written  for  an  oblate 
spheroid  are  applicable.  For  the  case  of  collapse  on  the  bottom,  the  cloud 
takes  the  shape  of  the  general  ellipsoid  and  a  frictional  force  between  the 
bottom  and  collapsing  cloud  is  included. 

Long-term  transport  diffusion 

11.  The  long-term  dispersion  phase  is  treated  in  one  of  two  ways.  When 
the  rate  of  horizontal  spreading  in  the  dynamic  collapse  phase  becomes  less 
than  an  estimated  rate  of  spreading  due  to  turbulent  diffusion,  the  collapse 
phase  is  terminated.  During  collapse,  solid  particles  can  settle  as  a  result 
of  their  fall  velocity.  As  these  particles  leave  the  main  body  of  material, 
they  are  stored  in  small  clouds  that  are  characterized  by  a  uniform  concentra¬ 
tion,  thickness,  and  position  in  the  water  column.  In  the  first  method  of 
handling  the  transport-diffusion  computations,  these  small  clouds  are  allowed 
to  settle  and  disperse  until  they  become  large  enough  to  be  inserted  into  the 
grid  positioned  in  the  horizontal  plane.  Once  small  clouds  are  inserted  at 
particular  grid  points,  these  points  then  have  a  concentration,  thickness,  and 
top  position  associated  with  them.  Figure  B2  illustrates  a  typical  concentra¬ 
tion  profile  at  a  grid  point.  Computations  on  the  grid  are  made  using  a  back¬ 
ward  convection  scheme  rather  than  attempting  a  numerical  solution  of  the 
governing  convection-diffusion  equation.  Tn  the  backward  convection  solution 
technique,  a  massless  particle  at  each  grid  point  at  the  present  level  is 
moved  backward  in  time  by  the  ambient  current  to  the  position  it  occupied  one 
time-step  before.  The  concentration  at  the  grid  point  it  presently  occupies 
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Figure  B2 .  Aspects  of  passive  diffusion 
(from  Brandsma  and  Divoky  1976) 

is  then  taken  as  a  five-point  average  of  points  surrounding  its  old  position. 

12.  Rather  than  making  computations  at  each  point  of  the  long-term,  grid 
at  each  time-step,  an  alternative  method  for  handling  the  transport-diffusion 
computations  has  been  incorporated  that  only  uses  the  horizontal  grid  for  out¬ 
put  purposes.  The  idea  for  this  method  was  obtained  from  work  by  Bra;,  isma  and 
Sauer  (1983),  on  the  development  of  a  drilling  mud  model.  Rather  tha^  insert¬ 
ing  the  mass  from  the  previously  discussed  small  clouds  into  the  horizontal 
grid,  each  small  cloud  is  assumed  to  have  a  Gaussian  distribution  given  by 
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where 

x,y,z  =  spatial  coordinates 

x  ,y  ,z  =  coordinates  of  cloud  centroid 
o  o  o 

o  ,o  >c  =  standard  deviations 
x  y  z  ^ 

m  =  total  mass  of  cloud  in  ft 

At  the  end  of  each  time-step,  each  cloud  is  advected  horizontally  by  the  input 
velocity  field.  The  new  position  of  the  cloud  centroid  is  determined  by 


x  =  x  +  u  ■  At 

o  o  , 

new  old 


z  =  z  +  w  •  At 

o  0  1  j 
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where 

u,w  =  local  ambient  velocities  in  ft/sec 
At  =  long-term  time-step  in  sec 

13.  In  addition  to  the  advection  or  transport  of  the  cloud,  the  cloud 
grows  both  horizontally  and  vertically  as  a  result  of  turbulent  diffusion. 
The  horizontal  diffusion  is  based  upon  the  commonly  assumed  4/3  power  law. 
Therefore,  the  diffusion  coefficient  is  given  as 

4/3 

Kx,z  =  \L 


where  is  an  input  dissipation  parameter  and  L  is  set  equal  to  four 

standard  deviations.  The  expression  for  the  horizontal  growth  of  a  cloud  then 
becomes 
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14.  Vertical  growth  is  similarly  achieved  by  employing  the  Fickian 
expression 
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where 

K  =  vertical  diffusion  coefficient 

y 

t  =  time  since  formation  of  the  cloud 
From  Equation  B5, 


and  thus, 
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-1/2  =  _y 


(B6) 


K 

a  =  a  +  - At  (B  7 ) 
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where  K  is  a  function  of  the  stratif ication  of  the  water  column.  The  maxi- 

y 

mum  value  of  is  input  as  a  model  coefficient  and  occurs  when  the  water 

density  is  uniform. 

15.  If  long-term  output  is  desired  at  the  end  of  a  particular  time- 
step,  the  concentration  of  each  solid  type  is  given  at  each  grid  point  by  sum¬ 
ming  the  contributions  from  individual  clouds  as 


where  N  is  the  number  of  small  clouds  of  a  particular  solid  type  and  y 
(the  vertical  position  at  which  output  is  desired)  is  specified  through  input 
data. 


B7 


16.  At  the  present  time,  the  effect  of  horizontal  solid  boundaries  has 
not  been  included.  Therefore  the  Gaussian  cloud  method  of  transport-diffusion 
computations  should  only  be  used  if  solid  boundaries  are  far  removed  from  the 
suspended  material.  However,  such  an  effect  due  to  the  bottom  and  the  water 
surface  has  been  included  in  the  vertical.  This  is  accomplished  through  re¬ 
flection  principles  by  assuming  that  identical  clouds  lie  above  the  water  sur¬ 
face  and  below  the  bottom. 

17.  In  addition  to  the  horizontal  advection  and  diffusion  of  material, 
settling  of  the  suspended  solids  also  occurs.  Therefore,  at  each  net  point 
the  amount  of  solid  material  deposited  on  the  bottom  and  a  corresponding 
thickness  are  also  determined.  A  basic  assumption  in  the  models  is  that  once 
material  is  deposited  on  the  bottom  it  remains  there,  i.e.,  neither  erosion 
nor  stable-bed  movement  of  material  is  allowed.  This  is  the  primary  theoreti¬ 
cal  limitation  of  the  models  that  restricts  their  usefulness  to  the  study  of 
the  short-term  fate  of  discharged  material. 

Model  Capabilities 


18.  The  computer  programs  enable  the  computation  of  the  physical  fate 
of  dredged  material  disposed  in  open  water.  The  following  discussion  de¬ 
scribes  particular  capabilities  or  special  features  of  the  codes. 

Ambient  environment 

19.  A  wide  range  of  ambient  conditions  are  allowed  in  model  computa¬ 
tions.  Conditions  (ranging  from  those  found  in  relatively  shallow  and  well- 
mixed  bays  and  estuaries  to  highly  stratified  two-layer  flow  fields,  found  in 
estuaries  where  salt  wedges  are  formed)  can  be  variable  from  one  long-term 
grid  cell  to  the  next.  The  only  restriction  on  bottom  topography  is  that 
associated  with  the  collapse  phase  which  was  discussed  in  paragraph  B5.  Any 
of  three  options  of  ambient  current  illustrated  in  Figure  B3  may  be  selected, 
with  the  simplest  case  being  the  time-invariant  profiles  shown  in  Figure  B3a 
for  a  constant  depth  disposal  site.  The  ambient  density  profile  is  input  as  a 
function  of  water  depth  at  the  deepest  point  in  the  disposal  site.  This  pro¬ 
file  may  vary  with  time  but  is  the  same  at  each  point  of  the  grid. 

Time-varying  fall  velocities 

20.  If  a  solid  fraction  is  specified  as  being  cohesive,  the  settling 
velocity  is  computed  as  a  function  of  the  suspended  sediment  concentration  of 
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a.  SIMPLE  ORTHOGONAL  VELOCITY  PROFILES  FOR  CONSTANT 
DEPTH.  APPLIED  EVERYWHERE  IN  FIELD 


b.  VERTICALLY  AVERAGED  VELOCITY  PROFILES  FOR  VARIABLE 
DEPTHS  WITH  EQUIVALENT  LOGARITHMIC  PROFILES 
SUPERIMPOSED. 


Figure  B3.  Illustration  of  the  various  velocity  profiles  available 
for  use  in  models  (from  Brandsma  and  Divokv  19761 


that  solid  type.  The  following  algorithm  is  currently  used 


I  =  0.0017  if  C  <  25  mg/i 

Vg  <  =  0.00713  C4/3/304.8  if  25  <  C  <  300  mg/?  (B9) 

(  =  0.047  if  C  >  300  mg/? 


where 

V  =  settling  velocity,  ft/sec 
s 

C  =  suspended  sediment  concentration,  mg/? 

Conservative  constituent  computations 

21.  The  models  allow  for  the  dredged  material  to  contain  a  conservative 
constituent  with  a  nonzero  background  concentration  of  that  constituent.  Com¬ 
puting  the  resultant  time-history  of  that  concentration  provides  information 
on  the  dilution  that  can  be  expected  over  a  period  of  time  at  the  disposal 
site . 

Output  available 

22.  Through  Input  data  the  user  specifies  the  amount  of  output  desired. 
Much  of  the  input  data  required,  e.g.,  the  water  depth  field,  are  immediately 
printed  after  being  read.  At  the  end  of  the  convective  descent  phase,  the 
location  of  the  cloud  centroid,  the  velocity  of  the  cloud  centroid,  the  radius 
of  the  hemispherical  cloud,  the  density  difference  between  the  cloud  and  the 
ambient  water,  the  conservative  constituent  concentration,  and  the  total  vol¬ 
ume  and  concentration  of  each  solid  fraction  are  provided  as  functions  of  time 
since  release  of  the  material. 

23.  At  the  conclusion  of  the  collapse  phase,  time-dependent  information 
concerning  the  size  of  the  collapsing  cloud,  its  density,  and  its  centroid 
location  and  velocity  as  well  as  conservative  constituent  and  solids  concen¬ 
trations  can  be  requested. 

24.  At  various  times,  as  requested  through  input  data,  output  concern¬ 
ing  suspended  sediment  concentrations  and  solids  deposited  on  the  bottom  can 
be  obtained  from  the  transport-diffusion  computations.  If  the  backward  con¬ 
vection  long-term  scheme  is  employed,  the  suspended  sediment  concentration  and 
the  location  of  its  "top  hat"  profile  (Figure  B2b)  are  provided  at  each  grid 
point  for  each  sediment  fraction.  However,  if  the  Gaussian  cloud  long-term 
scheme  is  selected,  only  concentrations  at  the  water  depths  requested  are 
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provided  at  each  long-term  grid  point.  In  both  cases,  the  volume  of  each 
sediment  fraction  that  has  been  deposited  in  each  grid  cell  is  provided.  At 
the  conclusion  of  the  simulation,  a  voids  ratio  specified  through  input  data 
is  used  to  compute  the  thickness  of  the  deposited  material. 

Assembly  of  Input  Data 

25.  Depending  upon  the  complexity  of  ambient  conditions  at  the  disposal 
site,  the  preparation  of  input  data  can  range  from  requiring  the  application 
of  a  three-dimensional  model  to  provide  stratified  velocity  fields  to  a  simple 
input  data  setup  of  perhaps  20  to  25  lines.  Input  data  can  be  grouped  into 
(a)  a  description  of  the  ambient  environment  at  the  disposal  site,  (b)  charac¬ 
terization  of  the  dredged  material,  (c)  data  describing  the  disposal  opera¬ 
tion,  and  (d)  model  coefficients. 

Disposal  site  data 

26.  The  first  task  to  be  accomplished  when  applying  the  model  is  that 

of  constructing  a  horizontal  grid  over  the  disposal  site.  The  number  of  grid 

points  should  be  kept  as  small  as  possible  but  large  enough  to  extend  the  grid 

beyond  the  area  of  interest  at  the  level  of  spatial  detail  desired.  Quite 
often  one  may  wish  to  change  the  horizontal  grid  after  a  few  preliminary  runs. 
Water  depths  and  the  horizontal  components  of  the  ambient  current  must  be  at 
each  grid  point.  Any  of  the  three  options  of  velocity  input  illustrated  in 
Figure  B3  may  be  selected,  with  the  simplest  case  being  velocities  at  a  con¬ 
stant  depth  disposal  site.  The  ambient  density  profile  at  the  deepest  point 
in  the  disposal  site  must  also  be  input.  This  profile  may  verv  with  time  but 

is  assumed  to  he  the  same  at  each  point  of  the  grid. 

C h nracterization  of  dredged  material 

27.  The  dredged  material  can  be  composed  of  up  to  12  solid  fractions,  a 
fluid  component,  and  a  conservative  chemical  constituent.  For  each  solid 
fraction,  its  concentration  by  volume,  density,  fal L  velocity,  voids  ratio, 
and  an  indicator  as  to  whether  or  not  the  fraction  is  cohesive  must  t>e  input. 
Proper  material  characterization  is  extremely  important  in  obtaining  rea'istic 
predictions  from  the  models.  If  a  conservative  chemical  constituent  is  to  be 
traced,  its  initial  concentration  and  a  background  concentration  must  he 
given.  Tn  addition,  the  hulk  density  and  aggregate  voids  ratio  of  the  dredged 
material  must  he  prescribed  along  with  its  liquid  limit. 
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Disposal  operations  data 

28.  Information  required  includes  the  position  of  the  barge  or  scow  on 
the  horizontal  grid,  the  volume  of  material  dumped,  and  the  loaded  and  un¬ 
loaded  draft  of  the  disposal  vessel. 

Model  coefficients 

29.  There  are  14  coefficients  in  DIFID.  Default  values  are  contained 
in  the  computer  code  that  reflect  the  model  developer's  best  guess.  However, 
the  user  may  input  other  values.  Computer  experimentation  such  as  that  pre¬ 
sented  by  Johnson  and  Holliday  (1978)  has  showrn  that  results  appear  to  be 
fairly  insensitive  to  man)  of  the  coefficients.  The  most  important  coeffi¬ 
cients  are  drag  coefficients  in  the  convective  descent  and  collapse  phases  as 
well  as  coefficients  governing  the  entrainment  of  ambient  water  into  the 
dredged  material  cloud. 
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